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The  annealing  behavior  of  defects  observed  in  electron  paramagnetic  resonance  (EPR)  and  photolumines¬ 
cence  (PL)  is  discussed.  We  consider  the  divacancy  (the  P6/P7  EPR  centers)  and  a  previously  unreported  EPR 
center  that  we  suggest  is  a  Vc-VSi-Vc  trivacancy  and  their  relationship  with  each  other  and  with  the  UD1-3 
series  of  PL  lines  near  1  eY.  We  observe  that  the  divacancy  and  the  UD2  PL  lines  annealing  behavior  is 
strongly  correlated  further  establishing  the  relationship  between  the  EPR  and  PL  centers.  We  present  a  detailed 
discussion  of  this  center  and  its  29Si  and  13C  hyperfme  spectra  which  supports  our  assignment  as  a  trivacancy. 
The  intensity  of  this  center  increases  with  annealing  temperature  as  the  divacancy  decreases  and  there  is  a 
sample-to- sample  correspondence  between  the  overall  intensities  of  the  two  centers.  However,  the  details  of 
their  annealing  suggest  a  more  complex  relationship  than  a  simple  one-to-one  transformation.  In  addition, 
while  the  UD1  PL  increases  with  annealing  temperature,  sample-to-sample  variations  indicate  it  is  not  related 
to  this  EPR  center. 
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As  wide  band  gap  semiconductor  technology  based  on  the 
Ill-nitrides  and  SiC  matures,  the  performance  of  high  power 
microwave  devices  based  on  these  materials  are  surpassing 
traditional  Si  technology.1  The  4 H  polytype  of  SiC  has 
emerged  as  an  especially  important  material  for  such  elec¬ 
tronics  because  of  its  availability,  high  electron  mobility, 
large  band  gap,  and  high  thermal  conductivity.  Point  and 
extended  lattice  defects  play  a  crucial  role  in  these  materials 
both  in  limiting  performance  and  in  compensating  shallow 
dopants  in  semi-insulating  (SI)  substrates.2,3  Structural  iden¬ 
tification  of  such  defects  is  important  to  understand  the  for¬ 
mation  of  these  defects  and  to  limit  or  optimize  their  concen¬ 
tration.  A  number  of  native  defects  such  as  the  isolated 
vacancies,  the  di vacancy,  and  the  carbon  vacancy-carbon  an¬ 
tisite  pair  have  been  identified  by  electron  paramagnetic 
resonance  (EPR).4  Although  the  identity  of  the  defect  or  de¬ 
fects  responsible  for  the  SI  behavior  or  limiting  device  per¬ 
formance  has  not  been  determined,  it  is  becoming  increas¬ 
ingly  apparent  that  more  than  one  defect  and  more  complex 
defects  may  be  crucial  to  understanding  the  compensation 
mechanism  in  SI  substrates  or  in  limiting  properties  such  as 
carrier  lifetime. 

There  have  been  a  number  of  experimental  and  theoretical 
studies  of  the  annealing  and  stability  of  fundamental  defects 
in  SiC.5-9  Much  of  the  experimental  work  has  been  on  irra¬ 
diated  material  where  a  primary  process  is  the  recombination 


of  vacancies  and  interstitials.  In  as-grown  material  there  are 
few  interstitials  and  so  such  processes  are  much  less  impor¬ 
tant.  In  such  material  the  C  vacancy  is  less  mobile  than  the  Si 
vacancy  with  the  former  being  stable  to  1600-1700  °C  (Ref. 
7)  while  we  find  the  Vsi  anneals  at  —1200  °C.  As  we  show 
later,  the  divacancy  is  also  quite  stable  to  —1500-1600  °C. 

In  this  work,  we  discuss  the  annealing  behavior  of  the 
di  vacancy  and  its  relationship  to  the  UD2  PL  lines.  In  addi¬ 
tion  we  focus  on  a  defect  (denoted  ANN1)  that  may  be  re¬ 
lated  to  the  divacancy  and  which  we  identify  as  a  trivacancy 
(Vc_^sr^c)  complex.  This  may  be  one  of  the  simplest  of  a 
family  of  more  complex  defects  that  play  a  role  in  the  SI 
character  of  these  substrates  and  the  starting  point  of  vacancy 
clusters  which  would  develop  into  more  extended  defects. 
EPR  work  has  also  been  reported  on  larger  less  well-defined 
clusters  of  vacancies.10 

I.  EXPERIMENT 

The  samples  used  in  this  work  are  standard  commercially 
available  physical  vapor  transport  (PVT)  or  high  temperature 
chemical  vapor  deposition  (HTCVD)  grown  SI  substrates 
with  B  and  N  levels  below  1016  cm-3  (Ref.  2).  Details  of 
sample  characteristics  before  annealing  are  given  in  Table  I. 
All  samples  were  annealed  for  30  min  in  an  argon  atmo¬ 
sphere  with  an  estimated  cool  down  time  of  —30  min,  which 


TABLE  I.  The  growth  techniques  and  initial  characteristics  of  the  sample  sets  used  in  this  study. 


Sample 

Growth 

Technique 

Initial  Characteristics 

Initial  Resistivity 

A 

PVT 

Strong  Si  Vac.,  Weak  C  Vac. 

1  X  106  fl  cm 

B 

PVT 

No  Si  Vac.,  Strong  C  Vac. 

5  X  1011  fl  cm 

C 

HTCVD 

No  Si  Vac.,  Strong  C  Vac.,  Strong  UD1  PL 

1  X  1011  fl  cm 

D 

PVT 

SI5(yc-Cs)  EPR 

5  X  1011  fl  cm 
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FIG.  1.  (Color  online)  A  representative  EPR  spectrum  illustrat¬ 
ing  the  P6/P7  centers  (the  divacancy)  in  an  unannealed  sample. 
Also  shown  are  the  isolated  vacancies. 

is  much  more  rapid  than  typical  cooling  rates  in  bulk  growth, 
possibly  quenching  in  metastable  defects.  All  EPR/PL 
samples  for  different  anneal  temperatures  from  a  given  wafer 
were  adjoining  pieces  cut  from  near  the  wafer  center.  At  the 
higher  temperatures,  the  annealing  resulted  in  visible  graphi- 
tization  of  the  surfaces.  The  thin  graphite  layer  was  removed 
by  annealing  in  an  oxygen  atmosphere  at  1000  °C  for  ap¬ 
proximately  one  h.  Electrical  measurements  on  sister  pieces 
annealed  at  the  same  temperatures  (but  a  different  time) 
showed  steady  decreases  in  the  room  temperature  resistivi¬ 
ties  with  annealing  from  ~106-1012  ohm -cm  to 
10-1000  ohm- cm.11 

The  EPR  spectra  were  obtained  using  a  conventional 
9.5  GHz  spectrometer  equipped  with  a  liquid  helium  flow 
system  for  temperature  control  with  optimal  data  obtained 
between  10  and  100  K.  The  EPR  studies  required  rapid  pas¬ 
sage  conditions12  and/or  optical  excitation.  Either  a  xenon 
lamp  and  various  filters  or  an  Ar+  laser  through  an  optical 
fiber  provided  the  required  optical  excitation.  Both  absorp¬ 
tion  in  quadrature  with  the  magnetic  field  modulation  and 
dispersion  measurements  were  used.  The  former  gives  a  de¬ 
rivative  line  shape  and  better  resolution  while  the  latter  gives 
a  stronger  undifferentiated  signal.  Relative  intensities  were 
obtained  by  double  or  single  integration  of  the  respective  line 
shapes.  Considerable  care  was  taken  to  maintain  the  same 
optical  intensity  from  run  to  run  and  results  were  quite  re¬ 
producible.  The  PL  measurements  on  the  same  pieces  used 
for  EPR  were  performed  at  1.6  K.  The  emission  was  excited 
by  the  351-nm  line  of  an  Ar+  laser,  analyzed  by  a  1/4-m 
double-grating  spectrometer,  and  detected  by  a  liquid 
nitrogen-cooled  Ge  photodiode. 

II.  RESULTS 

The  divacancy,  first  observed  by  Vainer  and  IP  in  in  irra¬ 
diated  6H-SiC  (Ref.  13)  is  commonly  seen  in  SI  substrates 


FIG.  2.  (Color  online)  The  photoluminescence  (PL)  spectrum 
illustrating  the  UD1,  UD2,  and  UD3  IR  PL  lines. 

and  has  been  associated  with  the  P6/P7  EPR  centers14-16  and 
these  in  turn,  have  been  linked  with  the  UD2  IR  PL  between 
1.1  and  1.16  eV.17  Magnusson  and  Janzen  have  recently  ar¬ 
gued  that  the  P6/P7  centers  are  associated  with  a  different  set 
of  PL  lines  between  1.0  and  1.06  eV.18  This  second  set  of 
lines  does  line  up  reasonably  well  with  the  MCD  (magnetic 
circular  dichroism)  detected  absorption  lines  which  have 
been  associated  with  the  P6/P7  EPR  lines  by  Lingner  et  al. 19 
However,  this  second  set  of  lines  is  only  seen  in  lightly 
ft-type  materials  while  UD2  and  P6/P7  are  common  to  SI 
material.  It  is  quite  possible  that  both  sets  of  lines  are  related 
to  the  di vacancy  but  in  somewhat  different  configurations. 
Here,  we  briefly  restate  our  original  arguments  and  provide 
further  data  in  support  of  the  coupling  of  UD2  and  P6/P7.  In 
Fig.  1  we  show  a  representative  P6/P7  EPR  spectrum.  The 
P6  lines  have  axial  C3V  symmetry  and  the  accepted  model  is 
a  Vc-VSi  divacancy  oriented  along  the  c  axis.14  The  P7  lines 
have  Clh  symmetry  and  are  due  to  a  divacancy  oriented 
along  a  basal  bond  direction.  In  Fig.  2  we  show  representa¬ 
tive  IR  PL  spectra  to  illustrate  the  UD1-3  lines  commonly 
seen  in  these  substrates.  Of  particular  interest  are  the  UD2 
lines  which  we  have  previously  associated  with  the  P6/P7 
EPR  centers.  Previously  we  have  shown  that  the  P6  centers 
have  sharp  zero  phonon  lines  (ZPL)  near  1.1  eV  in  photo- 
EPR  and  that  these  match  up  with  two  of  the  four  UD2  lines. 
The  P7  lines  do  not  have  sharp  ZPL  in  photo-EPR  but  their 
onset  does  match  up  reasonably  well  with  the  other  two  lines 
in  the  UD2  PL  spectrum.  In  this  work  we  discuss  the  anneal¬ 
ing  behavior  of  the  P6  centers  and  their  corresponding  UD2 
lines,  shown  for  two  samples  in  Fig.  3.  It  is  perhaps  not  too 
surprising  that  as  the  anneal  temperature  increases  both  the 
EPR  and  the  PL  spectral  intensities  decrease.  However,  note 
that  a  number  of  details  of  the  annealing  behavior  of  the  P6 
lines  are  reflected  in  that  of  the  UD2  PL  lines.  In  all  samples 
there  is  rapid  decrease  in  both  EPR  and  PL  between  1400 
and  1600  °C;  however,  in  sample  D  both  the  PL  and  EPR 
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Anneal  Temperature  (°C) 

1800  1700  1600  1500  1400 


FIG.  3.  (Color  online)  The  P6  EPR  (closed  symbols)  and  1.1  eV 
doublet  of  the  UD2  IR  PL  (open  symbols)  intensities  as  a  function 
of  anneal  temperature.  The  intensities  are  normalized  to  the  values 
of  sample  D  at  1400  °C. 


FIG.  4.  (Color  online)  The  EPR  spectrum  of  the  ANN1  center. 
The  up  arrows  indicate  the  A  =  15.4  MHz  29Si  hyperfine  satellites 
and  the  down  arrows  indicate  the  A =4.2  MHz  29Si  hyperfine 
satellites. 


level  off  and  decrease  only  very  slowly  between  1600  and 
1800  °C  while  in  sample  B  the  signals  drop  to  zero.  While  it 
is  possible  to  propose  a  number  of  models  to  account  for  this 
apparent  two  step  behavior,  the  fact  that  the  PL  mimics  the 
EPR  as  well  as  it  does  can  only  be  explained  by  both  being 
due  to  the  same  structural  defect.  Also,  we  observe  that  the 
slope  of  the  PL  and  EPR  intensities  vs  1/T  is  very  similar, 
especially  in  sample  B  (EA  ~  7  eV)  where  the  annealing  be¬ 
havior  is  rather  simple.  Finally,  we  remark  that  sample-to- 
sample  correlations  between  not  only  the  overall  intensities 
but  also  relative  peak  intensities  support  the  association  of 
the  UD2  PL  and  the  P6/P7  EPR.  In  summary,  the  annealing 
behavior,  the  optical  response,  and  sample-to-sample  varia¬ 
tions  all  support  the  identification  of  the  UD2  PL  as  arising 
from  the  P6/P7  EPR  centers  which  have  been  convincingly 
associated  with  the  divacancy. 

In  Fig.  4,  we  show  the  spectrum  of  a  defect  whose  par¬ 
entage  may  be  the  divacancy.  No  signal  is  seen  in  the  dark — 
this  may  be  because  this  defect  is  diamagnetic  in  the  ground 
state  but  can  be  optically  excited  into  a  paramagnetic  S  =  1 
triplet  state.  It  could  also  be  that  the  ground  state  is  the  triplet 
and  optical  excitation  serves  to  increase  the  population  of  the 
ms= 0  state  and  take  the  spin  system  out  of  thermal  equilib¬ 
rium  greatly  enhancing  the  intensity  of  the  EPR  lines.  It  has 
been  shown  that  this  is  true  with  other  defects  in  SiC. 14,20  In 
either  case,  the  phase  of  one  of  the  sets  of  lines  (the  low  field 
in  this  case)  are  inverted  because  the  center  is  being  optically 
pumped,  resulting  in  the  ms= 0  level  being  almost  com¬ 
pletely  filled  and  the  ms=±  1  levels  being  nearly  empty.21 
Hence,  the  low  field  lines  are  due  to  stimulated  emission 


while  the  high  field  lines  are  due  to  absorption,  similar  to 
that  shown  for  the  P6/P7  centers  in  Fig.  1. 

The  angular  dependence  of  this  spectrum  is  shown  in  Fig. 
5  for  rotations  about  the  [1120]  and  [0001]  directions.  For  an 
arbitrary  direction  of  the  magnetic  field,  there  are  three  mag¬ 
netically  inequivalent,  but  crystallographically  equivalent 
sites  of  the  defect  in  the  SiC  lattice.  When  the  external  mag¬ 
netic  field  is  restricted  to  a  {1120}  plane,  the  spectrum  con¬ 
sists  of  two  lines,  one  of  which  is  doubly  degenerate;  this 
degeneracy  is  removed  if  the  magnetic  field  is  out  of  the 
plane.  For  rotations  in  the  (0001)  plane,  the  spectrum  splits 
into  three  components.  Along  the  c  axis  all  the  sites  are  de¬ 
generate  and  only  one  pair  of  lines  is  observed.  The  points  in 
Fig.  5  are  the  experimental  data  and  the  solid  lines  are  gen¬ 
erated  using  best  fit  parameters  obtained  by  exact  diagonal- 
ization  of  the  matrix  form  of  the  following  spin-Hamiltonian 
(Ref.  22): 


H  =  g(3SB  +  D 


S2Z--S(S+  1) 
3 


+  e[s2x-sA  (l) 


where  j3  is  the  Bohr  magneton,  B  the  magnetic  field,  and  D 
and  E  are  the  axially  symmetric  component  and  the  asym¬ 
metry  parameter  of  the  crystal  field  interaction,  respectively. 
This  results  in  g  =  2.0056  (isotropic);  D  =  -279.8  MHz;  E 
=  5.6  MHz.  Clearly  the  principal  axis  of  the  crystal  field  in¬ 
teraction  is  along  a  (1100)  direction.  It  might  seem  unusual 
that  the  crystal  field  term  is  nearly  axial  even  though  the 
principal  axis  is  in  the  basal  plane.  However,  the  P7  defect 
(the  divacancy  oriented  along  a  basal  bond  direction)  has  a 
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Angle  (degrees) 

FIG.  5.  (Color  online)  The  orientation  dependence  of  the  ANN1 
center  for  rotations  about  a  (1120)  axis  on  the  left  and  about  the 
[0001]  axis  on  the  right.  The  positions  of  the  resonances  are  indi¬ 
cated  by  the  solid  circles  and  the  lines  are  a  fit  to  a  Hamiltonian 
containing  Zeeman  and  crystal  field  terms. 

range  of  \E/D\  values  from  0.2%  for  one  site  in  6H  SiC  (Ref. 
19)  to  20%  for  one  site  in  4 H  (Ref.  15)  and  so  the  ratio  of 
—2%  in  this  case  does  not  seem  unreasonable.  This  may 
indicate  that  the  defect  wave  function  is  highly  localized  and 
not  very  sensitive  to  the  surrounding  lattice.  Two  29Si  hyper- 
fine  values  were  also  measured,  A  =  15.4  MHz  and  4.2  MHz, 
each  due  to  between  two  and  four  29Si  nuclei.  The  orienta¬ 
tion  of  the  defect  suggests  a  pair  of  like  nearest-neighbor 
vacancies  while  the  intensities  of  the  Si  hyperfine  interac¬ 
tions  suggest  a  pair  of  C  vacancies.  The  29Si  hyperfine  split¬ 
tings  are  isotropic  to  within  experimental  accuracy  and  quite 
small  indicating  that  only  a  small  portion  (<10%)  of  the 
defect  wave  function  is  localized  on  those  atoms.  In  Fig.  6, 
we  show  a  simulation  of  the  hyperfine  structure  around  the 
high  field  line  assuming  the  4.2  MHz  interaction  is  due  to 
four  Si  neighbors  and  the  15.4  MHz  interaction  is  due  to  two 
Si  neighbors.  In  addition  to  the  primary  hyperfine  lines,  there 
are  weaker  lines,  as  indicated  by  the  arrows  in  Fig.  6  that  are 
due  to  two  29Si  atoms  around  the  center.  In  particular,  notice 
the  very  weak  lines  at  about  333.0  and  334.1  mT  in  the 
simulation.  These  are  due  to  two  29Si  atoms  on  the  two 
15.4  MHz  sites  and  are  too  weak  to  be  observed  above  other 
small  signals  in  the  data.  Other  possible  configurations  in¬ 
volving  three  or  four  atoms  on  these  sites  result  in  much 
stronger  lines  at  these  positions  which  would  be  visible 
above  the  background.  For  instance,  if  there  were  an  equal 
number  of  sites  with  both  the  4.2  and  15.4  MHz  hyperfine 
interaction  these  two  outer  lines  would  have  the  same  inten¬ 
sity  as  the  other  hyperfine  lines  due  to  two  29Si  atoms.  Such 
lines  would  clearly  be  above  the  background:  the  fact  that 
they  are  not  supports  the  model  with  only  two  atoms  with 
this  hyperfine  interaction.  In  Fig.  7  we  show  the  13C  hyper¬ 
fine  structure  that  we  believe  is  due  to  two  carbon  atoms 
which  are  nearest  neighbors,  along  with  the  two  C  vacancies, 

to  a  central  Si  vacancy.  For  B  ||  [1  f00]  the  hyperfine  interac¬ 
tion  is  84.2  MHz  and  the  intensity  of  these  two  lines  is 
—  1.6%  of  the  total,  consistent  with  two  C  neighbors.  By 
symmetry,  the  hyperfine  interaction  due  to  the  two  C  atoms 
should  be  degenerate  for  B||[1100]  and  should  broaden  and 
then  split  for  angles  away  from  this  orientation.  In  Fig.  7  we 


Magnetic  Field  (mT) 

FIG.  6.  (Color  online)  A  comparison  of  the  EPR  spectrum  of  the 
low  field  ANN1  line  and  a  simulation  assuming  two  Si  nearest 
neighbors  with  a  hyperfine  interaction  of  15.4  MHz  and  four  with  a 
hyperfine  interaction  of  4.2  MHz.  The  arrows  indicate  hyperfine 
lines  due  to  two  29Si  neighbors. 

observe  the  hyperfine  structure  does  split  for  small  rotations; 
although  the  resolution  is  not  as  good  as  one  would  like.  The 
individual  hyperfine  lines  appear  to  broaden  somewhat  and 
there  is  interference  from  other  weak  resonances  making  a 
full  rotation  study  infeasible.  The  principal  part  of  the  wave 
function  is  on  these  two  atoms.  Assuming  all  sp 3  wave  func¬ 
tions,  a  simple  atomic  orbital  calculation  yields:  10%  on 
each  of  the  two  C  atoms,  0.5%  on  each  of  four  Si  atoms,  and 
1.6%  on  each  of  two  Si  atoms.  If  the  carbon  atoms  relax 
away  from  the  Si  vacancy,  their  wave  functions  will  be  less 
s- like  and  the  calculated  concentrations  larger  and  the  data  in 
Fig.  7  suggests  significant  anisotropy  (e.g.,  for  only  a  10° 
rotation  the  two  pairs  give  interactions  of  74.1  MHz  and 
87.7  MHz  vs  82.7  MHz  and  85.9  MHz  expected  for  an  sp3 
orbital)  but  we  have  been  unable  to  obtain  sufficient  angular 
dependence  to  accurately  determine  the  full  13C  hyperfine 
tensor.  These  calculations  are  quite  simplistic  and  do  not  take 
into  account  effects  such  as  spin  polarization  which  can  be 
quite  significant.23 

We  originally  proposed  that  this  center  could  be  due  to 
two  carbon  vacancies  around  a  carbon  antisite:24  however, 
the  splitting  of  the  hyperfine  spectra  described  above  indi¬ 
cates  that  two  C  atoms  are  involved  and  we  are  not  able  to 
resolve  any  other  13C  hyperfine  lines  that  might  be  attribut¬ 
able  to  a  single  central  atom.  We  can  also  rule  out  a  Si  atom 
on  this  site,  since  the  hyperfine  lines  from  a  single  29Si  atom 
would  be  roughly  twice  as  intense  as  seen  for  the  two  13C 
atoms.  One  other  possibility  is  an  impurity  atom  with  no 
isotope  of  significant  abundance  (>1  at  %)  with  a  nuclear 
spin;  e.g.,  oxygen.  Lacking  evidence  to  the  contrary,  we 
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FIG.  7.  (Color  online)  The  13C  hyperfine  interaction  satellites 
around  the  low  field  EPR  line  of  the  ANN1.  The  angles  are  with 
respect  to  the  C  axis.  The  0=80°  spectrum  has  been  shifted  by 
about  1.15  mT  so  that  the  central  portions  of  the  spectrum  are  ap¬ 
proximately  aligned.  (Only  the  15.4  MHz  29Si  hyperfine  lines  are 
resolved,  the  4.2  MHz  lines  are  obscured  by  the  central  line.) 


adopt  a  provisional  model  of  the  defect  as  Vc^sr^c- 

There  are  two  possible  atomic  configurations  which 
could,  in  principle,  fit  the  data.  In  the  first,  the  two  carbon 
vacancies  lie  in  the  basal  plane  and  the  two  remaining  carbon 

atoms  lie  in  a  {1100}  plane.  In  the  second  configuration,  the 

roles  are  reversed  and  the  two  vacancies  lie  in  a  {1 100}  plane 
and  the  remaining  carbon  atoms  are  in  the  basal  plane.  The 
latter  is  preferred  since  the  spin  is  primarily  on  the  two  atoms 
in  the  basal  plane  aligned  along  a  (1100)  direction  which 
accounts  for  the  symmetry  of  the  crystal  field  term.  We  as¬ 
sume  that  the  observed  29Si  hyperfine  interactions  are  due  to 
the  six  Si  atoms  bonded  to  the  principal  pair  of  carbon  at¬ 
oms.  Naturally  other  assignments  are  quite  plausible  and  the¬ 
oretical  calculations  will  be  needed  to  sort  out  the  detailed 
assignments.  The  atomic  structure  of  these  configurations  is 
shown  in  the  “ball  and  stick”  pictures  in  Fig.  8  with  the 
principal  Si  and  C  atoms  labeled  (a),  (b),  and  (c).  Interest¬ 
ingly,  we  only  observe  a  signal  from  the  one  configuration, 
suggesting  that  the  configuration  with  the  carbon  vacancies 
in  a  {1100}  plane  and  the  other  two  C  atoms  in  the  basal 
plane  has  a  significantly  lower  formation  energy.  Further¬ 
more,  we  do  not  see  separate  signals  due  to  different  combi¬ 
nations  of  hexagonal  and  cubic  sites.  This  is  consistent  with 
the  wave  function  being  on  the  two  C  atoms  in  the  basal 
plane  with  weak  interactions  with  surrounding  planes. 

We  briefly  summarize  the  evidence  for  the  Vcr^sr^c 
model  for  this  center,  (i)  The  symmetry  of  the  crystal  field 
indicates  two  like  atom  sites,  (ii)  The  intensity  of  the  29Si 
hyperfine  interactions  points  to  a  pair  of  carbon  vacancies. 


FIG.  8.  (Color  online)  Models  showing  possible  configurations 
of  the  trivacancy.  The  model  on  top  with  the  two  carbon  vacancies 
in  a  {1 100}  plane  is  preferred  over  the  one  on  the  bottom  with  the 
two  vacancies  in  the  basal  plane  as  discussed  in  the  text.  The  larger 
balls  represent  Si  atoms  and  the  smaller  ones  C  atoms  with  light 
gray  balls  representing  the  vacancies.  The  (a)  Si  atoms  have  the 
4.2  MHz  hyperfine  interaction,  the  (b)  Si  atoms  have  the  15.4  MHz 
hyperfine  interaction,  and  the  two  C  atoms  marked  with  (c)  have  the 
principal  13C  hyperfine  interaction.  The  Si  site  marked  with  an  as¬ 
terisk  may  be  the  C  antisite  as  discussed  in  the  text. 

(iii)  The  13C  hyperfine  interaction  indicates  two  carbon  at¬ 
oms  are  involved  around  the  central  Si  site  between  the  two 
Vc’s.  (iv)  There  is  no  hyperfine  interaction  due  to  the  central 
site  which  is  most  probably  a  silicon  vacancy.  As  discussed 
in  the  previous  paragraph,  the  two  carbon  vacancies  are  in  a 
{1100}  plane,  while  the  other  two  carbon  atoms  around  the 
Vsi  lie  in  the  basal  plane. 

The  EPR  intensity  of  this  center  increases  significantly 
with  annealing  in  all  samples  as  well  as  showing  consider¬ 
able  dependence  on  the  starting  material,  as  shown  in  Fig.  9. 
We  note  that  samples  B  and  D  have  the  strongest  overall 
signals  and  also  have  the  largest  divacancy  signals  in  the 
unannealed  samples  suggesting  some  relationship  between 
the  di vacancy  and  ANN  1 .  The  Vc-VSi  -Vc  defect  could  result 
from  the  capture  of  a  C  vacancy  or  quite  possibly  a  Si  va¬ 
cancy  by  a  divacancy.  While  the  capture  of  a  C  vacancy  is 
straightforward,  the  mechanism  involving  the  Si  vacancy  re¬ 
quires  a  little  explanation.  If  a  Si  vacancy  moves  to  a  next- 
nearest-neighbor  position  to  the  Si  site  of  the  divacancy,  the 
C  atom  between  the  two  defects  could  move  into  the  Si 
vacancy  leaving  behind  a  C  vacancy  and  creating  a  C  anti¬ 
site.  This  is  similar  to  the  mechanism  proposed  by  Lingner 
and  co-workers  for  the  annealing  of  the  isolated  Si 
vacancy.19  The  intensity  of  the  divacancy  decreases  with  an- 
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FIG.  9.  (Color  online)  The  intensities  of  the  ANN1  EPR  signal 
as  a  function  of  anneal  temperature  for  the  four  samples  as 
indicated. 


nealing  as  that  of  the  Vc-^sr^c  defect  increases,  although 
there  is  not  a  one-to-one  correspondence.  If  it  were  a  simple 
conversion  of  the  divacancy  to  a  trivacancy,  the  sum  of  the 
two  intensities  would  be  roughly  a  constant.  As  we  see  illus¬ 
trated  for  sample  B  in  Fig.  10  this  is  hardly  the  case  and 
furthermore  we  see  that  the  concentration  of  the  divacancy 
decreases  sharply  before  that  of  ANN1  begins  to  increase 
and  has  dropped  below  the  detection  limit  while  the  ANN1 
signal  continues  to  increase.  The  ANN1  signal  decreases  af¬ 
ter  higher  temperature,  presumably  due  to  conversion  to  an¬ 
other  defect  which  we  are  not  able  to  detect  in  EPR.  Similar 
results  are  observed  for  the  other  samples.  Previous  reports 
have  shown  the  concentration  of  carbon  vacancies  does  not 
change  appreciably  with  annealing  up  to  1700  °C  while  the 
Si  vacancy  becomes  mobile  at  —1200  °C.  While  we  do  not 
detect  the  Si  vacancy  after  these  anneal  temperatures,  only 
the  negatively  charged  state  has  been  identified;  although, 
the  positive  and  neutral  states  are  thought  to  be  in  the  gap 
and  paramagnetic.  One  might  speculate  that  the  process  is 
something  like: 


Vc-Vsi+ysi^Vc-Vsi-C-Vsi^Vc-ysi-Vc-Csi 

(2) 


in  which  the  intermediate  structure,  Vc-Vsi-C-Vsi,  is  diamag¬ 
netic  or,  at  least,  not  observed  in  our  EPR  measurements. 
The  C  antisite  would  be  reasonably  removed  from  the  prin¬ 
cipal  atoms  in  the  defect  wave  function  as  shown  in  Fig.  8 


Anneal  Temperature(°C) 

FIG.  10.  (Color  online)  The  normalized  intensities  of  the  diva¬ 
cancy  and  ANN1  signals  and  their  sum  for  sample  B.  The  divacancy 
is  normalized  to  its  value  at  TA=  1400  °C  and  ANN1  to  its  value  at 
TA=1800  °C.  Similar  results  are  observed  for  other  samples. 


and  should  relax  away  from  the  vacancy.  Therefore,  it  should 
have  little  effect  on  the  EPR  signal.  This  would  account  for 
the  minimum  in  the  sum  of  the  divacancy  and  ANN1  EPR 
signals  around  TA~  1550  °C  as  shown  in  Fig.  10. 

The  IR  PL  was  monitored  as  a  function  of  anneal  tem¬ 
perature.  The  UD1  PL  signal  as  shown  in  Fig.  11  increases 
with  anneal  temperature  while  the  UD2  and  UD3  PL  signals 
decrease.  A  cursory  look  at  Figs.  9  and  11  might  suggest  a 
relationship  between  UD1  and  the  ANN1  EPR  center.  How¬ 
ever,  while  both  the  UD1  PL  line  and  the  ANN1  EPR  center 
increase  more  or  less  monotonically  with  anneal  tempera¬ 
ture,  UD1  increases  by  a  factor  of  about  8  vs  a  factor  of 
about  3  for  the  EPR.  Furthermore,  there  is  little  correlation 
between  samples  (e.g.,  sample  C  has  the  strongest  PL  signal 
and  one  of  the  weaker  EPR  signals,  while  samples  B  and  D 
have  strong  EPR  and  weak  UD1  PL.)  It  seems  more  likely 
that  UD1  and  UD3  are  inversely  related,  although  we  do  not 
have  a  good  model  for  either.  This  is  illustrated  in  Fig.  12  in 
which  we  plot  the  sum  of  the  intensities  of  the  UD1  and  UD3 
PL  signals  normalized  to  their  values  after  1800  °C  and 
1400  °C  anneals,  respectively.  The  three  PVT  samples  (A,  B , 
and  D)  give  a  sum  which  has  a  standard  deviation  of  9-13% 
about  the  mean.  The  HTCVD  sample  (C)  follows  a  some¬ 
what  different  behavior,  perhaps  because  it  is  grown  at  a  very 
different  temperature.  We  see  in  Fig.  11  that  the  intensity  of 
UD3  in  sample  C  increases  until  it  reaches  a  maximum  after 
the  1550  °C  anneal.  (UD3  decreases  monotonically  with  an¬ 
neal  temperature  in  the  other  three  samples.)  If  we  replot  Fig. 
12  in  which  the  UD3  signal  has  been  normalized  to  its  value 
after  the  1550  °C  anneal,  the  sum  is  then  relatively  constant 
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FIG.  11.  (Color  online)  The  intensities  of  the  UD1  PL  as  a 
function  of  anneal  temperature  for  the  four  samples  as  indicated. 
Note  the  correspondence  between  symbols  and  samples  is  the  same 
as  in  Fig.  9  (e.g.,  sample  A  is  squares  and  so  forth.) 

between  1550  °C  and  1800  °C  anneals.  Both  UD1  and  UD3 
change  by  about  an  order  of  magnitude  with  annealing  while 
their  sum  changes  by  only  about  +10%,  suggesting  they  are 
related.  Magnus  son  and  co-workers18  have  argued  that  these 
centers  may  be  due  to  impurities  and  their  inverse  relation¬ 
ship  might  further  argue  that  they  are  due  to  the  same  impu¬ 
rity;  e.g.,  in  different  charge  states. 

III.  SUMMARY 

High  temperature  anneals  were  used  to  study  the  evolu¬ 
tion  of  native  defects  in  semi-insulating,  ultrahigh  purity 
4H  SiC.  We  present  additional  evidence  to  support  the  as¬ 
signment  of  the  UD2  PL  lines  to  the  P6/P7  EPR  lines,  which 
are  due  to  the  divacancy.  The  UD2  PL  has  very  similar  an¬ 
nealing  behavior  to  that  of  the  divacancy,  even  mimicking  a 
two  step  behavior  seen  in  some  samples.  We  also  see  that 
sample-to- sample  variations  are  consistent  with  the  relation¬ 
ship  and  previous  work  had  shown  that  the  zero  phonon 
absorption  lines  of  the  P6  center  observed  in  photo-EPR 
aligned  with  the  emission  line  observed  in  PL.  We  observe  a 
defect  that  we  preliminarily  identify  as  the  Vc_^sr^c  triva- 
cancy  based  on  the  symmetry  of  the  crystal  field  term  and  the 
29Si  and  13C  hyperfine  interactions.  The  largest  hyperfine  in¬ 
teraction  is  with  the  two  C  atoms  next  to  the  Si  vacancy. 
Based  on  the  symmetry  of  the  crystal  field  splitting  we  be¬ 
lieve  these  C  atoms  are  in  a  basal  plane  while  the  three 
vacancies  lie  in  a  {1100}  plane.  The  intensity  of  this  center 
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FIG.  12.  (Color  online)  The  sum  of  the  UD1  PL  normalized  to 
its  value  after  an  1800  °C  anneal  and  the  UD3  PL  normalized  to  its 
value  after  a  1400  °C  anneal  as  a  function  of  anneal  temperature. 
Data  for  samples  A,  B,  and  D  have  been  offset  for  clarity.  The  arrow 
indicates  +10%  deviation  from  the  mean.  Also  included  is  the  data 
in  which  the  UD3  PL  for  sample  C  (the  HTCVD  sample)  were 
normalized  to  its  value  after  a  1550  °C  anneal  (open  triangles). 
Note  the  correspondence  between  symbols  and  samples  is  the  same 
as  in  previous  figures  (e.g.,  sample  A  is  squares  and  so  forth.) 

increases  with  annealing  temperature  for  all  sample  sets  and 
is  strongest  in  samples  that  had  the  strongest  initial  diva¬ 
cancy  signals.  It  is  proposed  that  the  di vacancy  anneals  by 
capturing  an  isolated  vacancy  (either  C  or  Si).  However,  a 
simple  capture  mechanism  would  imply  a  one-to-one  corre¬ 
spondence  between  the  divacancy  and  the  trivacancy.  Since 
this  is  not  observed,  we  suggest  an  intermediate  diamagnetic 
species  which  would  account  for  the  minimum  in  the  com¬ 
bined  EPR  signal  around  1550  °C.  The  UD1  defect  observed 
in  IR  PL  also  increases  with  annealing;  however,  neither  a 
sample-to-sample  comparison  nor  the  details  of  the  anneal¬ 
ing  behavior  of  the  EPR  center  and  UD1  support  a  relation¬ 
ship  between  UD1  and  ANN1.  For  instance,  the  samples 
with  the  strongest  ANN  1  EPR  signal  (B  and  D)  have  weak 
UD1  PL  while  the  sample  with  the  largest  PL  (C)  has  a  weak 
EPR  signal.  Our  results  are  consistent  with  a  relationship 
between  the  UD1  and  UD3  PL  signals  which  indicate  that 
both  are  due  to  the  same  impurity  or  lattice  defect. 
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